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1^ , Abstract 



We calculate sin^ 9 (Mz) in the MSSM in terms of oem, Gp, 

rrit, tan/3 and SUSY mass parameters with the same accuracy as 

the present calculations of sin^ 6{Mz) in the SM. We compare the 

5^ \ results with the standard leading logarithmic approximation used for 

SUSY threshold corrections and find important differences in the case 
of light sparticles. We give approximate formulae connecting coupling 
constants in the SM and in the MSSM and comment on process depen- 
dence of such formulae. The obtained values of the MSSM couplings 
ai{Mz) are used to investigate gauge coupling unification in the min- 
imal SUSY SU{5) model. Our non-logarithmic corrections lower 
the predicted value of the Higgs triplet mass. The interplay between 
the supersymmetric and GUT thresholds in achieving unification for 
the coupling constants in the range of the experimentally acceptable 
values is quantified. 
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1. INTRODUCTION. 

The succesful gauge coupling unification in the Minimal Supersymmetric 
Standard Model (MSSM) [|I| has spurred recently a lot of interest in the SUSY 
GUT scenario. The standard approach to the coupling unification is based on 
the fact that physics at a given scale can be described in terms of an effective 
renormalizable theory, with the effects of the larger scales decoupled. For 
SUSY GUTs, one naturally considers the three energy regimes described by 
the Standard Model (SM), the MSSM and the full GUT. The first two models 
are effective renormalizable theories which approximate the full theory up to 
higher dimension operators suppressed by the scales of new physics. 

The standard approach to the gauge coupling unification consists of the 
following steps: 

i) Measurements at energy E < Mz are used to extract the couplings^ 
ai{Mz) i = y,2,3 using the Dimensional Reduction regularization (DR) [Q 
and the Modified Minimal Subtraction Scheme {MS 0) in the framework 
of the SM. In addition to a^i^Mz) and aEM^Mz) it is convenient to use 
sin^ 6w{Mz). It can be calculated e.g. in terms of Gp = 1.16639(2) x 10~^ 
GeV-^ a'ilj = 1/137.0359895(61), Mz = 91.1888(44) GeV §, rrit 
and M^o (the top quark and SM Higgs boson masses respectively). The 
first three parameters are known with very high accuracy and the latter 
two enter through radiative corrections. In practice, one takes aEM^Mz), 
sin^^^^(M^) and a^i^Mz) determined in the ordinary Dimensional Reg- 
ularization (DIMR) and MS and subsequently converts them to the DR 
regularization using the prescription |^ 

^' (1) 



ai{Mz) af^^'^^{Mz) 127r 

with Cy = 0, C2 = 2, C3 = 3. The Weinberg angle can also be extracted 
(or crosschecked) from other observables: e.g. Gp can be replaced by 
sin^^^ measured at LEP and SLD [^. However, at present the most 
precise determination of the couplings in the SM follows from the measured 
values of Gp, Mz and apM- 

a) The coupling constants are evolved to higher scales by means of the 2- 
loop RGE of the SM, up to some scale Msusy, and then using the RGE of 
the MSSM. The threshold corrections to the running of ai{Q) due to the 
splitting in sparticle masses are included at the 1-loop level and result in the 
contribution ^, |^, ^ 



1 1 TTli' 

^MSSM(M V - ^SM(M T = 47r A, (M5[/5y) = E ^^^^ log — 



^We call Uy the coupling related to the ordinary hypercharge reserving the symbol ai 
for the rescaled coupling appropriate for discussion of unification; ai = (5/3)ay. 



where the coefficients Abik are the contrbutions to the (3i function of 
the couphng ai from the sparticle k. It has to be stressed, that the 
effects of the VEV contribution to sparticle masses as well as the effects of 
the left-right mixing between sfermions cannot be treated consistently in this 
way and must be neglected. Thus the masses ruk appearing in the formula 
(0) are the lagrangian SU{2) x U{1) invariant mass parameters and not the 
physical masses. 

The one loop formula for a^^^^^{Q) at some scale Q > max{mk): 

^ - ^ + 26riog.^"^^^ 



aMSSM^Q) af^'{Mz) ' ^ M, 



z 



,MSSM 1 _ „„ Q 



+ 47rA,(M5c/5y) + 26f"^^'Mog— ^ (3) 

^^■isUSY 



can be rewritten as 



l l , oaf5SAfi„„ Q 



aP^*^(Q) af^^^^(Mz) ^ " M, 



+ 26f ^^^^ log -^ (4) 



with 



1 



af^^^iMz) af^iMz) 



+ 477 A, (Mz) (5) 



having the obvious interpretation of the couplings aj(M^) extracted at 
the scale Mz directly in the MSSM. This is possible, because the threshold 
corrections are given in eq.(2) in the form in which the 1-loop logarithms are 
resummed. Therefore the dependence on Msusy in eq.(3) drops out and 
one can take Msusy = Mz- One should stress that this is true only in the 
leading logarithmic approximation for the threshold corrections. In general 
the results would depend on the scale Msusy and this dependence is always 
a next order correction. 

Hi) The MSSM equations are used for the running up to some scale Mqut 
where the unification is achieved (or not) without or with the GUT scale 
threshold effects which are taken into account also in the leading logarithmic 
approximation. In the latter case the final result depends on the specific 
GUT. In addition it depends on our criterion of an "acceptable" GUT spec- 
trum. 

The standard approach is justified in the presence of well separate scales 
{Mz -C Msusy ^ Mgut)- However, it often happens that the scales 
of "new physics" are not large enough to neglect its non-logarithmic con- 
tributions. This may well be the case with the Standard Model and its 
super symmetric extension: it is quite plausible that the superpartner masses 
vary from 0{Mz) (some of them can even be lighter than Z° boson) to, 
say, (9(1 TeV). In this case the (renormalizable) SM is not a consistent 



approximation at the scale Mz to the full MSSM and higher dimension 
(non-renormalizable) operators cannot be neglected in step i), i.e. in the 
procedure of extracting the SM couplings from the data. 

An equivalent, more straightforward and more convenient approach is 
to work at the scale Mz directly with the MSSM. This point has been 
emphasized by Faraggi and Grinstein |TD|. This means that the couplings 



at Mz extracted from experimental data are the MSSM (and not the SM) 
couplings. All SUSY threshold effects are already included at this step and 
the RG running is based on 2-loop MSSM equations from Mz up to the 
GUT scale. From this point of view formula (|^) is a formula for a^^^^{Mz) 
calculated e.g. in terms of Gp, Oi^M, Mz, rrit and other SUSY parameters 
with only 1-loop logarithmic SUSY contributions retained, and resummed to 
all orders. We call it Leading Logarithmic Threshold (LLT) Approximation. 
The main part of this paper is the calculation in Sec. 2 of the sin^ 9]y (Mz) 
(and also, as a byproduct, of Myy) in terms of «^f/, Gp, Mz, rrit, tan/5 
and other SUSY parameters in the MS Scheme directly in the MSSM 
with the same accuracy as the present calculations of siii^ 9 w{Mz) in the 
SM. In Sec. 3 we compare the results with the standard LLT approximation 
and find important differences in the case of light sparticles. Also in Sec. 3 



we comment on the Farragi-Grinstein (FG) |T0| approach and its ambiguity 
related to the process dependence. In an approximation analogous to theirs 
we give the formulae connecting the coupling constants in the SM and in the 
MSSM, if calculated in terms of Gp, Mz and a^f/ . In Sec. 4 we use the 
obtained values of the MSSM couplings ai{Mz) to investigate the gauge 
coupling unification in the minimal SUSY SU{5) model. In particular, we 
find that our non-logarithmic corrections lower the predicted values of the 
Higgs triplet mass. Also, we study in some detail the interplay between the 
super symmetric and GUT thresholds in achieving unification for the coupling 
constants in the range of the experimentally acceptable values. 

2. CALCULATION OF sin^ OwiMz) IN THE MSSM. 

The calculation of sin^ 9w{Mz) proceeds through the calculation of the 
Fermi constant measured in the /x — >■ eUgUfj^ decay. We follow closely the 
analogous calculation in the SM whose details are analyzed in [|ll[]. The 
equation for s^ = sin^ 9]y(Mz) reads (for the running coupling constants 
in the MSSM taken at the scale Q = Mz we use the abbreviations: a = 
apAiiMz) and dj = ai{Mz) for i=y,2,3; analogous couplings in the SM 
will be distinguished always by the superscript "SM".): 

V2GpM^l-Afw{a,Mz,Mw,mt,...) ^^ 



where 



Afw 



n 



WW 



(0) 



M^ 



+ 






+ 5 



VB 



(7) 



and M]y (Mz) is the physical mass of the W^ (Z°) boson. The 11^-^^2(5^) 's 
are the vector boson self energies calculated in DR and renormalized by the 
Modified Minimal Subtraction 



MS- 



and Z'-y self energies in DR can be found in [T^) and 6vb = ^v^ '^^vW^ 



the expressions for WW, ZZ, 77 

>VB — OvB 

is the contribution of "non-oblique" vertex, box and external fermion wave 
function renormalization factors. In the MS scheme, the SM contribution 
(5y^ is given by |TT| : 



"VB 



a 

Atts' 



6 + 



0„2 



10), 2 

logCiy 



(8) 



with c^ = 1 — Syy used as an abbreviation for M^^/M^- It has the same 
form in the DIMR and DR schemes. The Sy^^^ has exactly the same 
form as in the on-shell calculation (with the obvious replacement of sw, cw 
and a^ff by s, c and a) and can be taken over from the Appendix A 
of ref. 
in the 



m 



the counterterms to the vertices and neutrino self energies have 
MS calculation the interpretation of the electron and muon wave 



function renormalization factors. 

One way to organize the calculation with Gp, Mz and a^f/ taken as 
the input parameters is [|TT| to express the physical mass of the W^ boson 
in terms of Mz'- 



M'w 



c'pMJ 



(9) 



with 



1- 



Uww{M^) 



1 




M^ 


p 


1- 


iizz{Ml) 5tlzz{Ml) 
Ml Ml 



+ 5p 



(10) 



where 6Ilzz{Ml) is the contribution of the Z — 7 
between the physical and bare Z^ boson masses: 

2 



5f[zz{Ml) 



{tiz.jMl 

Ml-n, 



,,(M|) 



mixing to the relation 



(11) 



and Sp = 5fP'~^^ + Sp^^'^'-^^ stands for the leading irreducible 2-loop 
corrections to p of order 0(03^^^?) (QCD) and C(G|.m^) (HIGGS). 
6p^^^ in the SM and in the limit of a heavy top quark reads [IT 



Sp 



QCD 



2a2.{mt) 




(12) 



where rrit is the top quark pole mass. This expression is expected to account 
for the dominant part of QCD corrections also in the MSSM since corrections 
coming from gluons attached to squark loops, which themselves contribute 



-HIGGS 



was 



a small piece of p, should be very small. In the SM, Sp 
calculated in ref [15| in the limit M^^o = 0. Recently this calculation has 
been extended to M^o 7^ by Barbieri et al. [|T^]. In the absence of explicit 
calculation of ^pHiGGS -^^ ^^^q MSSM we use the following interpolating 
formula: 



5p 



HIGGS 



-Nr 






vr^ 



cos a 
sin/3 



i(2) 



M, 



hO 



nit 



+ 



I sma 



\ sin /3 ^ 



i(2) 



M, 



HO 



nit 



(13) 



where p^'^\x) is given in |T6[ (see also '^^), M^o and Mj^o are masses of 
the lighter and heavier MSSM scalar Higgs bosons, tan/3 = f2/t'i and a 
is the mixing angle between scalar Higgs bosons. 

Finally, the running coupling constant a = a{Mz) is obtained from the 
formula [p.?}] : 



a 



a 



Ad 



(14) 



with a = a^f/ and 



7a , Mw 8a , rrit 
Ad = 0.0684 ± 0.0009 + — log — ^ log tt^ 

Att ^ Mz 97r ^ Mz 



Ana A, 



(15) 



where 



Stt^A, 



24 mc^ I Mh± 

> - log H — log 

^3 ^ Mz 3 ^ Mz 



+ > - log 



+ > - log 

^9 ^ Mz 



1=1 



J^4:^ Mu- 



contains logarithms of the physical masses of all charged supersymmetric 
particles. Notice, that as compared to the DIMR (see e.g. ^7^), the quan- 
tity Ad in DR does not contain any constant term coming from integrating 
out W^ and Z° bosons f\ 

The term 0.0684 ±0.0009 accounts for the contribution of light fermions 
and contains a nonperturbative part. It is extracted from the experimental 
data on e+e~ — > 7* ^ hadrons |jl^. The error 0.0009 is the most important 
theoretical uncertainity in the determination of s^ from eq. (E). 



^This is the content of the theorem S stating that in the DR scheme there are no 
threshold corrections when gauge bosons are integrated out. 



The fine structure constant a^f/ is obtained from tlie Tlionison scat- 
tering in pure QED, wliicli is tlie effective theory at the scale of the electron 
mass. At this scale, all higher dimensional operators which remain after 
integrating out all heavier particles are totally negligible and so are the cor- 
rections to eq.([T^). 

Eq.(^ is solved iteratively by calculating at each step of iteration the 
physical W boson mass from eq. @ and a from eq. ( |14|) . 

In the limit of very heavy sparticles SUSY corrections to self energies 
of the vector bosons are dominated by large logarithms of the soft SUSY 
breaking masses. In this limit our approach should give the same results 
as given in eq.(^ in the LLT approximation, where the logarithms read 
explicitly: 

3 ^ 1 3 



stt^a^ = -yiog— ^ + -yiog— ^ 

1 J^ mn- 4J^, rriu. 1 J^, mn- 

+ - y bg — ^ + - y bg —^ + — y bg — ^ 

9 ^ ^ Mz 9^, ^ Mz 18^, ^ Mz 

2, |/i| 1 Mao .^^. 



8^ ^2 = 7T 2. bg -— + - }_^ bg 



^U ^Mz 2^ -Mz 

+ - log — log — log — ^ f 18 ) 

3 ^ Mz Q Mz 3 Mz ^ ' 

( mQ^,mL., ... are the soft SUSY breaking mass parameters). From eq.(|^) 
one easily gets: 

^ ^ +47r(A^ + A2) (19) 



a a^^^ 



.2 ., 1 + 47raf ^A, 



where quantities with the superscript SM are determined by equations 
^) in the SM i.e. with the contributions of SUSY particles removed. 



Analogous formula with af^^A2 replaced by Oy^^Ay holds for c^. Since 
in the limit of heavy sparticles the physical masses are dominated by large 
soft SUSY breaking masses we have Ag = A„ + A2 + 0( ..2^ ) and it is 

easy to see that indeed a calculated from eq. (0) approaches a determined 
from eq.(|l9D. Next, it is important to notice that in the limit considered, by 



using the asymptotic form of the SUSY contributions to self energies Q: 



flzz{q')-tlzz{0) ~ A7i^q^{c'A2 + s^Ay), 

tiz^iq^) ~ Any^q^ (c^A^-s^Ay), (21) 



n^^(g2) ~ 47rd g2 (Aa + Aj,). 

eq.(P) can be rewritten in the following form: 

£2_^2 1 + 47rdf ^A, 1 (00) 

''^' 1 + 4™^^^ (A, + A2) 1 - AfsM47rdf ^ A2 ^ ^ 

which is the same as eq.(|l9]) (up to higher order terms). Thus, we reproduce 
correctly the leading logarithmic threshold corrections in the limit of heavy 
sparticles. 



With regard to eqs.(^ and (10) we see that they connect with each other 
two observables Mw and Mz- Therefore, according to the Appelquist- 
Carazzone decoupling theorem, heavy SUSY particles have no effect on this 
relation; one can explicitely check this cancellation in the limit of heavy 
sparticles with the help of eqs.(21) and (p^). 

Before presenting our results for sin^ 9{Mz) (see Sec. 3) we comment 
on various consistency checks of our calculation. The higher order correction 
in eq. (^ coming from neglected (non-top exchange) 2-loop diagrams can be 
estimated as in the SM [^ to be of order 6s'^ < 3 x 10^^; the uncertainty 



due to the difference between eqs.(^) and (0) turns out to be < 1 x 10 '^ for 
MsusY < 1-5 TeV. The W^ mass calculated in the MSSM in WS scheme 
(eqs.(|) and (p!OD) and in the on-shell scheme^ |13[ agree to better than C(10) 
MeV for a wide range of sparticle masses. The W^ mass calculated in the 
MSSM (by using eqs.(|^) and (PUJ)) with sparticle masses greater than 1.5 
TeV agrees to better than 10 MeV with the W^ mass calculated in the 
SM (in the MS scheme or in the on-shell sheme) with the appropriate SM 
Higgs boson mass M^o. 



■^Subtraction of n(0) for W^ and Z° self energies is required because of the 
chargino/neutralino contrubution which unhke scalars do not decouple in that limit from 
renormalized self-energies at q'^ = 0; however, Ilww{0)/M^r — Ilzz{0)/M^ vanishes 
in the limit of heavy charginos and neutralinos as it should. 

"'in the present version, our code used in |13|| has been improved to incorporate 
gpHiGGS g^g jj-^ eq.(p^) and the leading top dependent corrections to the Higgs boson 
masses using the effective potential approach [n9| as described in EG]. 



3. COMPARISON WITH THE LEADING LOGARITHMIC 
AND FARAGGI-GRINSTEIN APPROXIMATIONS. 



We begin the presentation of our results with a brief summary on s|j^^ 
{siji^ 9w{Mz) in the Standard Model) calculated with our code in terms of 



Gp, M. 



a 



a 



OS 

EMi 



rrit and M^o. A comparison with earlier results 
ll|, [l^ is a useful check of our calculation and remembering several numerical 



values provides a convenient reference frame for comparison with the MSSM 



results. 
Mw: 



In Table 1 we collect a sample of values for s|jy/ (in DIMR ) and 









Table 1 






rrit 


160 


160 


170 


170 


180 


180 


M^o 


60 


130 


60 


130 


60 


130 




0.23155 


0.23195 


0.23125 


0.23164 


0.23093 


0.23132 


Mw 


80.344 


80.301 


80.406 


80.362 


80.471 


80.426 



We reproduce the results of ref. |Tl| within ±0.00008 accuracy in s|^.^ for 
a wide range of rrit and M^o values. Our input parameters (with which the 
Table 1 has been obtained) are: Mz = 91.1888 GeV, Gp = 1.166739 x 10~^ 
GeV-2, a = a% = 1/137.0369895 and a^iMz) = 0.118. The running 



electromagnetic coupling d" 



in the MS scheme reads: 



a 



1 8 ^ ruf 

3P = 127.87+- log ^ 



(23) 



and is calculated from eq.(|14|) by taking 



, Of 7a , Mw 8a , rrii 

Ad = 0.0684 + — log — ^ log TT^ 

Qtt 477 ^ M^ 97r ^ Mz 



(24) 



appropriate for the calculation in DIMR (the conversion to DR can always be 
performed with the help of eq.(|l|)). The uncertainty in the value of a^'^ re- 
lated to the hadronic contribution to Ad ((5(Ad) = ±0.0009 |]18|) is a 
source of uncertainities: Ss"^ = ±0.0003 and 6Mw = ±17 MeV. 

Similarily, our values of the physical mass Mw obtained in the Standard 
Model from eqs.(^) and ([l0| ) agree very well with the previous calculation [11, 
|T^ (up to a few MeV). As mentioned earlier, the present calculation in the 
MS scheme also agrees with the calculation in the on shell renormalization 
scheme: for 90 < rrit < 240 GeV and 60 < M^o < 1000 GeV the 
difference is less than 5 MeV and for rut < 210 GeV and M^o < 500 GeV 
the agreement is typically within 2 MeV. 

^Provided as in ref. we take Mz=91.17 GeV, SpQ'=^" = 0, 
and adjust the hadronic vacuum polarization to their value. 



Sp 



HIGGS _ 



= 27r^ - 19, 



9 



For later use it is very convenient to have simple interpolating formulae 
for s|^ and Mw in the SM. For the central values of the input parameters 
we find that the following formulae: 

s|^ = 0.23166 + 5.4 x IQ-^h - 2.4 x IQ-^h^ 

- 3.03 X lO^^t - 8.4 X lO^V (25) 

{Mw)sM = 80.347 - 6.4 x 10~^h + 2.5 x IQ-^h"^ 

+ 6.2 X lO^^t + 1.25 X 10"¥ , (26) 

where h = M^o — 100 and t = m^ — 165 (both masses in GeVs), reproduce 
the results for s|^^ and M\y in the SM with the accuracy 0.00001 and 2 MeV 
respectively (for the range 130 < rrit < 200 GeV and 60 < M^o < 140 GeV 
which is relevant for comparison with MSSM). These formulae are therefore 
used as input in eq.(^ in our further study whenever we refer to the values 
of s|jv/ and My^/ in the Standard Model. 

Turning now to our results for s^ in the MSSM, there are two points to 
be discussed first: the departure from the SM result (with the same values of 
the input parameters Gp, Mz, a^M, rrit and M^o) as a function of the 
SUSY particle masses and comparison with the LLT approximation P, ^]. 
Both are illustrated in Fig.l a-d where we present the results of our complete 
calculation and of the LLT approximation for a scan over the sparticle masses 
in a wide range of values given in the Appendix. In our scan we include 
only those values of the parameters which are consistent with experimental 
constraints on the sparticle masses, the mass of the lighter MSSM scalar 
Higgs boson |^ and give Mw within 2a of the measured value [^]. In 



Fig.l, plots for rrit = 180 GeV contain smaller number of points because 
most of the cases with light sparticles are eliminated by the mentioned above 



experimental cuts |[T3| , |20[| . All points for nit = 160 (180) GeV in Fig.l 
correspond to the values s^^ = 0.23155 - 0.23195 (0.23093 - 0.23132), 
with the ranges reflecting the (weak) dependence of the s|^,^ on M^o in 
the range 60 — 130 GeV (see Table 1) f\ and show the dependence of the 
running couplings in the MSSM on SUSY thresholds. 

Comparing the complete calculation with the LLT approximation we see 
that the difference in s^ is up to (9(0.002) for light supersymmetric spec- 
trum and therefore very relevant for the discussion of unification of couplings. 
There is also a nonnegligible dependence of a on the contribution from the 
electroweak symmetry breaking to the charged sparticle masses (see Figs. 2- 

4). 

The origin of the differences between the full and approximate calculations 
can be understood from formulae quoted in the Appendix. Qualitatively, the 

^This dependence is neglected in the frequently used fit s|j^^ = 0.2326 — 1.03 x 
lO-^K - 138^) 0. 

10 



largest deviations from the LLT approximation come from the contribution 
of sfermion and, to the smaller extent, chargino/neutralino sectors. This is 
further illustrated in Figs. 2-4 where we show the dependence of our results 
on the chosen sparticle masses keeping other SUSY particles heavy (C(l 
TeV)). 

Finally we discuss the approximate method proposed by Faraggi and 
Grinstein [^ of translating the values of the couplings in the SM into the 
values of the couplings in the MSSM, without performing the complete cal- 
culation in the MSSM, as in Sec. 2. 

In this approach one is instructed to choose a set of observables Oi 
which can be calculated in the SM and MSSM as 



O. = F,{gt 



SM\ 



(27) 



and 



O,: 



fp ( MSSM 



(28) 



respectively, where e„ stand for SUSY parameters not present in the 



SM. Writing next g^^^^^ 



^F,{g^ 



MSSM 



9l'' + 59k 



H9, 



MSSM 



H9, 



MSSM\ 



+ 



Cn) and equating the RHS of the above two equations one ar- 



rives (assuming that 6gk are small enough) at the set of equations 



^2^9 



QgSM 



+ AF,igt',e. 



(29) 



which allow for the determination of the MSSM couplings provided the SM 
couplings g^^^ have been determined from the chosen set of observables Oi. 
This prescription can be easily applied to our set of observables: Q^em, 
Gp and Mz leading to the following approximate formulae (which can 
be also derived directly from eqs.(6-10) without any reference to the FG 
method): 



5gl 1 

2 ^2 ^2 
92,SM ^SM ~ ^SM 




nzz(M|)\"^"''" 

Ml ) 


^9l 1 

2 "2 ^2 
9y,SM ^SM ~ '^SM 




- / , .9 X \ SUSYl 

^zz{Ml)\ 

Ml I \ 


or, equivalently 







(30) 



5a 



a 



SM 



-47rd^*^Ae 



(31) 
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Ob _ Lgj^j 



"2 "2 "2 



:SM^ , ^I^WwiO) UzziM, 



^ . SUSY- 



Ana''''' A, + Tie 



M^ Ml 



(32) 



Similarily, one can easily find the formula for correction to tlie W^ mass 
{5 Ml, = iM^)MssM - {M^)sm)- 



Tie 



{M^)sM c|m \sImJ V M^ Ml 




In the above formulae, (...)'^^'^^ stands for superparticle contributions to the 
vector boson self energies and includes also the difference between the Higgs 
sector contributions to these self energies in MSSM and SM. The formulae for 
tlwyviO)^^^"^ and Ilzz{Ml)^'-'^^ are collected for the reader's convenience 
in the Appendix. 

Using the interpolating formulae (|25| , ^6D for s|^j and {M]y)sM "we can 



check that equations (|3T| , p2D give an excellent approximation to our complete 
results for the couplings. This is illustrated in Figs. 2-4. The approximation 
for {Mw)mssm given by eq.(^) is typically within 20 MeV as compared 
with the full calculation in the MS scheme. 

One should stress that the FG approximation is consistent only if the 
SM couplings are determined from the same set of observables Oi as used 
to calculate the couplings in the MSSM: Qk'^^^^ = 9k^^ + ^Qk- This is 
because the corrections 5gk are process dependent and correlated with 
the values of g^^^ extracted from the same set of observables: since we 
are concerned with the situation when the SM (without higher dimension 
non-renormalizable operators) is not a sufficiently accurate effective theory 
at Mz scale, the values of the couplings f^f*^ extracted from different 
experiments (a different set of observables Oi) are also process dependent 
and g^^{Oi) — gl^^{Oj) can be of the order of 5gk itself. It is therefore 
inconsistent to apply the FG approach to s\^ obtained from an overall fit to 
the data in the SM; such a fit gives meaningless numbers if non-logarithmic 
SUSY corrections are important. It is also inconsistent to apply FG method 
to correct the SM couplings obtained from a set Oi by using the equations 
for 5gk derived by referring to a set Oj. This is shown in Figs. 2-4 where 
the results obtained from the two variants of the FG formulae [|10| derived by 



reffering to the scattering processes and applied to the SM couplings obtained 
from our set of observables M^, Gp and a^ff are marked by dotted 
lines. 

It is easy to check, using formulae (21) that all the approximations give 
the same asymptotic behaviuor of the corrections to a and s^ when the 
SUSY particles are heavy. Nevertheless, as is clear from Figs. 2-4, for light 
SUSY particles they give very different results from the complete calculation. 

4. UNIFICATION OF COUPLINGS. 
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Our calculation of the running Weinberg angle in the MSSM can be used 
to investigate unification of the gauge couplings, with the supersymmetric 
threshold effects taken exactly into account at the 1-loop level. The main un- 
certainty in our value of sin^ 9 w{Mz) (for fixed m^, Mf^o and superpartner 
masses) follows from the uncertainty in a which gives (5s^ = ±3 x 10~^ (the- 
oretical uncertainity of the calculation like higher order corrections etc. is 
smaller, O{lxl0~^)). This can be compared with the precision ±6 x 10~^ 
of sin^ 9 w{Mz) in the SM obtained from the global fits to the LEP and 
the low energy data |2^, |^. Note, however, that as explained earlier, from 



the point of view of the MSSM the value of the sin^ 9 w{Mz) in the SM 
obtained from global fits is meaningless as soon as we want to go beyond the 
LLT approximation. To increase the statistical significance of the value of s^ 
in the MSSM, i.e. to obtain it from a global fit to more observables than just 
Gp, Mz and a^f/ one has to perform similar calculation to ours for other 
observables and to have an independent fit for each chosen set of values for 
rrit and M^o and for superpartner masses. In particular it is important to 
remember about the correlation between the value of sin^ 9 w{Mz) and the 
values of the top quark and Higgs boson masses which is already present in 
the Standard Model results (see e.g. ref. [0 and eq. (H)). 



Given the high accuracy of our approximate formulae (|31|j3^ , |33|) which 
translate s|j^,f into s^ in the MSSM (for any chosen superpartner spectrum) 
and hoping for similar accuracy of the formula connecting af ^^ with a^ 
in the MSSM 0]: 



,\SUSY 
&3^ "^ \ Ml 



i^_^,(IW^. (3„ 



where H^^'^^ is the SUSY contribution to the gluon self energy (provided 
we apply this formula to af ^ extracted from processes where the gluon 
four- momentum q^ is of order of M|), it is convenient to use the values 
of s|^ and af^ as a label for the set of values of P and a^ in 
the MSSM which correspond to the same values of s|^j and af^ in the 
approximation (|3T| , |32| , |33| , |3^ . The values used in the following are as in Table 
1: 

mt = 160 GeV sIa/ = 0.23155 - 0.23195 
mt = 180 GeV s|m = 0.23093 - 0.23132 

For the strong coupling constant af*^ we shall be using values 0.115, 
0.125 and 0.13, for the purpose of illustration of various effects. Thus, 
in this paper, by "experimental data" we mean the above values of s|j,j, 
df*^ and the corresponding sets of values of s^, ds in the MSSM (now 
dependent on SUSY parameters). Exactly analogous remarks apply to a'^'^ 
and a. The "experimental data" for the MSSM will be compared with the 
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relation between s^, a and ds predicted after imposing the unification 
conditions in the MSSM. Note that since we work directly in the MSSM, 
the predictions which follow from unification can be obtained, as a relation 
between s^, 0:3 and a without any reference to the "experimental data". 
The final point we have to specify before the actual discussion of unifi- 
cation is the choice of the supersymmetric spectrum. The very general scan 
over the sparticle masses performed in Section 3 is not very illuminating from 
the point of view of unification: by stretching the parameters and correlat- 
ing them in some specific way one has very little constraint on unification 
scenarios. It is much more interesting to study unification of the gauge cou- 
plings under the additional assumption that the supersymmetric spectrum is 
obtained from the requirement of radiative electroweak symmetry breaking. 



In this paper we use the spectra of ref. p^] obtained with universal boundary 



conditions at the GUT scale and with one - loop corrections to the effective 
potential included. The spectra of ref. [^] correspond to a complete scan- 
ning over the available parameter space with cut-off on the physical squark 
masses Mg < 2 TeV. The requirement of radiative breaking results in very 
strong correlation between masses of different sparticles. 

From the point of view of the gauge couplings unification a useful param- 
eter which characterizes the spectrum is Tsusy defined so that in the LLT 
approximation (neglecting GUT threshold corrections) the SUSY threshold 
correction to the value of af ^ predicted from the gauge coupling unification 
reads 0, §: 

A(l/d3) = --log^^^ (35) 

Assuming that soft SUSY breaking mass terms are the same for all genera- 
tions we get (see the Appendix for notation): 
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"--'"'l^J [4) l#J [^) '''^ 

One of the consequences of the correlations in the sparticle spectra ob- 
tained in the model with radiative electroweak breaking is that, generally, the 
effective parameter Tsusy < 300 GeV for Mg < 2 TeV, and it is strongly 
correlated with the supersymmetric parameter fi (only in the fixed point 
regime, Tsusy extends to 600 GeV for Mg < 2 TeV and this is due to 
large values of /i, which are in this case necessary for proper radiative 
electroweak symmetry breaking) [^, |25|. This is shown in Fig.5 for several 



values of rrit [] and tan/3. It is important to realize that, in the context 
^The mass rrit is the pole mass which is related to the running mass as follows: 

mt = mtimt) (l + ^^ + 11-4 (^)') 
14 



of radiative electroweak symmetry breaking with universal boundary condi- 
tions at the GUT scale, the often referred to value of Tsusy = 1 TeV would 
correspond to very heavy squark spectra. 

We are now equipped to study unification of the gauge couplings. First, 
we address the question of unification with no GUT threshold and higher 
dimension operator GUT corrections. 

This requirement, as discussed before, determines a relation between s^, 
ds and a (all three are MSSM couplings). Fortunately the "experimental 
value" of d changes only weakly for the spectra of ref. [|4|: {<ymin{max))~^ 



= 131.8 (129.0) for heavy (light) charged sparticles. Therefore we present 
the unification prediction (obtained from 2-loop running with the MSSM 
RGBs from Mz to the GUT scale) as a narrow band in the s^ — ds plane 
obtained for a in its "experimental range". This band can be compared 
with "experimental data" for s^ and a^ obtained with the use of the 
chosen SUSY spectra, bearing in mind that {amin{max))~^ is more appro- 
priate for heavy (light) spectra. Labelling the "experimental data" by the 
corresponding Standard Model values of s|j\^ and df^ we see in Fig. 
6 a clear trend: unification with no GUT threshold corrections occurs only 
for sufficiently large df^ and sufficiently heavy supersymmetric spectrum 
(the actual numbers depend on nit and tan/3 values). The dependence on 
rrit has two sources: firstly it enters into the 1-loop correction to s^ and 
also to ds and d as already in the SM (see the values (pS])). Secondly 
the features of the SUSY spectrum depend on nit and tan/? because of 
the tight constraints imposed by universality of the soft terms and by ra- 
diative symmetry breaking. In particular, as discussed above and seen in 
Fig. 5 for the same cut - off for Mg the maximal value of Tsusy is 
nit and tan/5 dependent. In Fig. 7 we plot the predicted df*^ in the 
approximation ( PT| , |52|J5^ (from unification with no GUT thresholds) as a 
function of Tsusy for several values of nit and tan/5 and with s'sm 
and d^M in their experimental ranges specified earlier. The uncertainty 
6s'^ = ±3 X 10~^ (correlated with uncertainty in d, eq. (p!5|)) translates itself 
into Sa^^ = =f1-5 x 10~^ in this plot. The dependence of s|;y^ on nit 
(eq. (p5D) is of course refiected in the predicted value of df*^. Changing 
nit from 160 GeV to 160 ± 20 GeV (and leaving the sparticle spectrum 
unaffected) results in a shift dt&l^^ = ±(2 - 2.5) x 10"^ . 

In the LLT approximation and neglecting the dependence of s'^j^^j on 
the Higgs boson mass, the predicted df^ is a function of Tsusy only. 
The scatter plot obtained in the LLT approximation reflects the dependence 
of the input s|^j on the Higgs boson mass which is varied from 60 to 
130 GeV. In the full calculation there is stronger departure from a universal 



In this paper, as in ref. [Q, we neglect other contributions to the pole mass. 
®The coupling constant usually quoted in the literature which we call df *^ throughout 
the paper is actually dg in the theory with the top quark decoupled. 
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dependence of af ^^ on Tsusy which is due to the non-logarithmic threshold 
corrections. Going beyond the LLT approximation results in larger values 
of df^^ for small Tsusy- Moreover, in the whole range of Tsusy 

reachable with our spectra, the corrections to the LLT approximation are 
non-negligible. This happens because the spectra always contain some light 
sparticle(s) which is responsible for the difference. In the absence of the GUT 
threshold corrections the unification scale can be unambigously defined as 
the point of crossover of the couplings. It depends only weakly on the SUSY 
spectra and varies between 1.3 x 10^^ and 3 x 10^^. 

Next we study unification with GUT thresholds included. We restrict 
ourselves to the minimal supersymmetric SU{5) with three different mass 
scales: M^s, M^ and My corresponding to the Higgs triplet, Higgs 
bosons originating from the 24 representation and heavy vector bosons 



(leptoquarks) respectively. We have then ||2l], |26 




aM) a.W)'2.V*'°^^-"°"l^-'°"^, 
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for any scale Q close to My^.z/s i-e. a scale at which the higher order 
(logarithmic) threshold corrections to the RHS can be neglected. It follows 
from eqs. (^) that M^^, and a combination of M^ and My are uniquely 
related to the values of the coupling constants: 
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(38) 

and ai{Q) are related to dj = ai{Mz) by the RGE in the MSSM. 

The curves of fixed Mh^ in the (s^, as) plane are given by eq. (^) 
(for fixed cxem)- They are shown in Fig. 6 for Mh^^ = 5 x 10^^, 2 x 
10l^ 1 X 10^^ (for aEM = 1/131.8). Thus, each point in the (s^ag) 
plane corresponds to fixed values of Mh^ and of the product MyM-£ 
(modulo the uncertainity in cxem)- Running with 2-loop equations from 
Mz to Q such that LHS of the second of eqs. (|38|) vanishes we obtain 
(M^Ms)^/^ (= Q) and Mh^. With GUT threshold corrections, unification 
of couplings can be always achieved at the expense of large enough splitting 
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between Mho, and the combination (My-Ms)^/^. There are, however, certain 
natural requirements which constrain these three masses and their sphtting. 
First of all, Mh^-, Ms and My should be smaller than the Planck mass. 
Other requirements follow from the superpotential which breaks SU{'b) and 
the condition that its couplings remain perturbative up to the Planck scale. 
This gives [|21|: 



Mh; < 2 My, Ms < 1.8 Mv (39) 

Another requirement is more technical: the discussion of unification based 
on the 1-loop threshold corrections eqs. (|3^) is meaningful only if there exist 
a common scale Q such that all logarithms \og{MH-i/Q), \og{MY,/Q) and 
\og{Mv/Q)-, are small enough to make higher order corrections (next to 
leading logarithms) negligible. Those arguments still leave a lot of freedom 
for GUT corrections. It is interesting to know how much, say, the predicted 
af^ is modified by the splitting {M^My)^/^ /Mh^ » 1. Using the 
formulae ( |38D with 1-loop running to Q = Mz we obtain 

which yields 5gut&I^'^ = -0.01 (-0.02) for {M^M^y/yMH, = 30 (1000) 
(which is not strongly modified by 2-loop effects) respectively. We demon- 
strate the effect of GUT scale threshold correction in Fig. 6 where we show 
the bounds in the (5^,0:3) plane in the MSSM which are consistent with 
unification with 

max(MH^, Ms, My) , , 

. ;,/^' ,, ' ,/( < 30 or 1000 41 

min{MH,, Ms, My) ^ ^ 

respectively. The only (weak) dependence of those bounds on the super- 
symmetric spectrum is through the value of oiem- The asymmetry of 
the bounds with respects to the unification curve with no GUT thresh- 
old corrections, seen in Fig. 6, is due to condition (|39|) . Below this curve 
My = Ms = 30 (1000) Mh:, and above M^g = 2My = 30 (1000) Ms which 
results in a much smaller splitting between M^^ and (M^Ms)^'^^. 

As stressed in ref. ^^, the uncertainty in the unification prediction for 
da due to the unknown GUT thresholds makes direct tests of the unification 
idea impossible, even if we knew the low energy superparticle spectrum and 
no matter how precise the experimental data for the couplings at Mz are. 
Only the reversed programme is realizable: the more precise the experimental 
data for the couplings are the more accurate predictions we can obtain for 
Mffg and for the splitting between Mh^ and (My-Ms)^''^. Comparing the 
bounds with our "experimental data" it is clear that for light SUSY spectrum 
GUT threshold corrections, corresponding to large splitting My(s)/M/^3 ~ 
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30 are needed for unification with df*^ = 0.125. If one takes af*^ = 0.115, 
close to its experimental lower bound, the same GUT splitting is needed 
even for Tsusy = 300 GeV and it has to reach My(Y.~^/MH^ ~ 1000 for 
light spectra. For the same value of af^, the non- logarithmic corrections 
calculated in this paper tend to increase GUT corrections needed to achieve 
unification. 

It is also very interesting to show the predictions for Mh^ which follow 
from our "experimental data" points. The Higgs triplet mass is directly 



related to the proton lifetime ^\- There have been claims in the literature 
that even Mh-j, ~ 10^^ GeV can be reconciled with the gauge coupling 
unification and at the same time it is large enough to make the proton life 
time consistent with experimental limits |^ also for large values of tan/3 



[ pl[] . Our predictions for Mh^ are plotted in Fig. 8 as a function of Tsusy 
for df*^ = 0.125. The value of Mh^ is rather sensitive to changes of the 
MSSM couplings at Mz- For example, a change 5s^ = 3 x 10~^ increases 
Mh^ by a factor of ~1.5 whereas a change 5a^ = —2 x 10^'^ decreases it 
by a factor ~2.2. 

We stress again that those results depend on the SUSY spectrum which is 
taken from the minimal supergravity model (radiative electroweak symmetry 
breaking with universal boundary conditions at the GUT scale) and with cut- 
off on the physical squark masses, Mg < 2 TeV. We obtain an upper bound 



on Mff^ which is lower than the one in ref. |21]. For the same value of 



/i = 1 TeV, df*^ = 0.125, rrit = 160 GeV and including the uncertainty 
(5s^ = ±3 X 10^*^ we obtain M^^^ = 3 x 10^^ . The discrepancy with 
the upper bound given in ref. ||21| is due to a different s'gj^^ value used in 



that paper. Our bounds on M^^ are in agreement with those given in ref. 



2^. We also observe that our non-logarithmic SUSY threshold corrections 
decrease the predicted value of Mh^. 

Finally we would like to discuss the impact on the unification of gauge 
couplings of the condition that quark Yukawa couplings remain in the per- 
turbative regime. In particular we impose 

F/ < 4vr (42) 

This condition is relevant only for a heavy top quark and small tan/5 (i.e. 
close to the quasi IR-fixed point 0, §]). As it is clear from the structure of 
the RGE for the Yt, for given values of rrit and tan/5, the condition (^) 
gives a lower bound on a^^^^^^ f\ Superimposed on the unification plots 
in Fig. 6, it favours light superpartner spectrum: for heavy spectrum the 

^The bound depends (weakly) on sin dw because a change in sin 6\y induces a 
change in the final point of the renormalization evolution where Ytop grows very quickly. 
In this discussion we neglect the contribution of SUSY and GUT threshold corrections 
to the Yukawa couplings and in this sense our conclusions here should be taken as only 
qualitative. 



values of a^^^^^^ corresponding to resonable values of df*^ (say, below 
0.13) are too small to be consistent with perturbative top quark Yukawa 
coupling. It is interesting to observe that, for small tan/?, large mt, 
and af*^ < 0.125(0.115), perturbativity of Yt puts strong upper bounds 
on TsusY- TsusY < 200( 600) GeV for m^ = 160 (180) GeV and 
tan/? = 1.25 (2.0) respectively. 

5. CONCLUSIONS. 

In this paper we have calculated the running Weinberg angle at Mz 
directly in the MSSM in terms of Gp, Mz, dEM, f^t and supersymmet- 
ric parameters. The accuracy of this calculation is similar to the analogous 
calculation in the Standard Model. Using the language of supersymmetric 
threshold corrections to the Standard Model our calculation is equivalent 
to a complete calculation of those corrections at the one - loop level. A 
detailed comparison with LLT approximation is presented. Also highly ac- 
curate approximate formulae are given which translate the SM couplings into 
the MSSM couplings. 

Those results are subsequently used to study the gauge coupling unifi- 
cation in the minimal supersymmetric SU{5) model, with the superpart- 
ner spectra taken from the minimal supergravity model (i.e. with radiative 
electroweak symmetry breaking and universal boundary conditions for soft 
supersymmetry breaking parameters at the unification scale). 

For unification with no GUT thresholds, the value of the strong coupling 
constant is calculated as a function of supersymmetric spectra. The non- 
logarithmic threshold corrections increase the predicted value of «f*^ for 
light spectra. For the spectra obtained in the minimal supergravity model 
with radiative electroweak breaking, with cut-off Mq < 2 TeV, one gener- 
ically obtains af^ above the experimental range af*^ = 0.118 ± 0.007 
(with the exception of the low tan/3 values and rrit close to its quasi- 
IR fixed point). In the presence of GUT thresholds, a direct precise test of 



the unification idea in the minimal SU(5) is no longer possible [26]. Only 
the reversed programme is realizable: with precise information on the cou- 
plings at Mz we can calculate the Higgs triplet mass and the combination 
(MyMj^Y^^ as a function of the supersymmetric spectrum. For the same 
values of the couplings at Mz, the non-logarithmic corrections calculated 
in this paper decrease the value of Mh-j and increase the splitting between 
Mh^ and (M^M^)^/^ needed to achieve unification. Information from the 
proton life time (sensitive to Mh^^) and from b — t Yukawa coupling unifi- 
cation (sensitive to the splitting) gives additional constraints on the model. 
It will be interesting to use the results of the present paper for improving the 
precision of those constraints. 
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APPENDIX. 

In this Appendix we make explicit our convention and notation for spar- 
ticle masses. In order to avoid confusion with signs, often present in the 
literature, we will be explicit here. We follow closely the convention (but not 
the notation) of pO|] . 

From the superpotential of the model 

W = eijYiHlUEt, + eijYdH{QW+ + eijYuHiQ^UR + ei,^iH\Hi (43) 

(where Ya are the Yukawa couplings and ei2 = — e2i = — 1) and the relevant 
soft SUSY breaking part of the lagrangian 

Csoft = - m^J^iT - ^hJH^I^ + mleij{H\Hi + c.c.) 

- mi|Lp-m||E+p-m||Q|2-m^|f/^p-m^|D+P (44) 
+ e,,YiAiH{UE+ + ei,YdA,H{QWl + e^^YiA^Hi^Q^U^ 

where Hi, H2 and L, Q are the SU{2) doublets and adding contributions 
of the D-terms one gets the mass matrices for up and down squarks as well as 
charged sleptons (because we neglect the intergenerational mixing we write 
them for one generation only; therefore U (u) represents generically up-type 
squarks (quarks) and the same is understood for dow-type ones and sleptons 
(leptons)): 



with 



Mfj 




Mk 
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ml + ml + ^t{Ml-AM^) 


ML 
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ml + ml- \t{Ml - Ml) 


Mlr^n 
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—■mu{Au + ficot/3) 


Mk 
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ml + ml~ h:{Ml + 2M^) 


Mk 
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ml + ml + U{Ml - Ml) 


Mkn 
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—nid{Ad + /itan/3) 


Mk 
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ml + ml- h{Ml - 2Ml) 


Mk 


= 


ml + ml + t{Ml - Ml) 


Mkn 


= 


—rrieiAe + fi tan j3) 



(45) 



(46) 



(47) 



(48) 
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Here t = (tan^/? — l)/(tan^/3 + 1) and mu,d,e stand for ordinary fermion 
masses. These matrices are diagonalized by the appropriate rotations. The 
sneutrinos have masses given by: 

Ml=ml + ml + ^M| (49) 

Part of the lagrangian relevant for the chargino/neutrahno sector reads: 

C = ^Mg,rr + lMg,xX-f^e,,h\hi + tV2g2HlT-hr (50) 
+ iV2g2H*T-h2r - -^gyH^hx + -^gyH^h^x + h.c. 

where tl)"" (a=l,2,3) and x are the SU{2) and U{1) gauginos respectively, 
hi and /i2 are the two higgsino doublets and T"" are the SU{2) generators. 
Definig v^V^^ = "^^ =F ^"^^ we get in the basis {—i^^,hl,—iip~,h'l) = 
(x+, x") the chargino mass matrix: 

with 

y_f Mg2 V2Mwsmp\ , . 

^~\V2Mw cos/3 /i J ^^^^ 

Diagonalization requires two unitary mixing matrices Z+ and Z_ which 
rotate x^ fields to mass eigenstates A^: A^ = Z}^x'^- The two physical 
chargino masses are denoted as mc\, i = 1,2. 

The neutralino mass matrix in the basis x^ = {~'^X) —i'4'^-, h\, h"^ has 
the well known form: 



M 



N 



( Mgi -Mzs COS (3 Mzs sin [3 \ 

Mg2 Mzc COS p -Mzc cos p 

-Mzs cos P Mzc cos P — /i 

V Mzs sin P -Mzccosp -fi / 



(53) 



and is diagonalized by rotation A*^ = Z^^x°. Physical neutralino masses 
are denoted as rrifyf., i = 1,...4. Four component chargino (Dirac) and 
neutralino (Majorana) fields are then built as 

Finally, gluino, which does not mix with anything, has a mass denoted 
as Mg. 
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We do not describe here the Higgs sector referring the reader to refs. 
[|1^, |12[. We recall only that at the tree level, the Higgs boson masses and 
couplings can be conveniently parametrized by the ratio of the vacuum ex- 
pectation values of the two Higgs doublets tan/? = f2/f i and the CP— odd 
Higgs boson mass My^o. 

In terms of the parameters defined above, the general scan used in Sec. 3 
can be described as follows. 50 < M40 < 500 GeV, 50 < mc , rn^. < 1500 
GeV, 120 < Mb^ < 1000 GeV, OAMb^ < Mt^ < 2.5Mt^, o' < At < Mb^, 

and An = for 



Qn = Mb, 



with Ai = (for all the 



OAMb, < Mi, < 0.8Mb, keeping Mq^ = M( 
the first two generations as well as Mej^ = My 
three generations). 

Next, we display the formulae for {Iiww{(f)Y^^^ and (Yizz{q^)Y^^^ 
which allow easy use of our approximation given by eqs. (|25|p6| , |3T|j3^ , |33| ) . As 
explained in Sec. 3, besides the genuine sparticle contributions, {J1ww{q^)Y^^^ 
and {JlzziQ"^))^^^^ include also the difference of contributions coming from 
the Higgs sectors of MSSM and SM and therefore depend also on the SM 
Higgs boson mass M^o. Obviously, calculating the corrections to s^ one 
has to take the same M^o in formulae below and in the fit (p5|). It is 
convenient to choose M^o = M/^o. 

Let us first consider W^ self energy. The difference of the MSSM and 
SM Higgs boson contributions gives: 



2\\h 



^Ti{Iiww{q')) 



-,Ml, 



bo{q\Mw,M^o) 
cIMq^ Mw, Mho) - slM<l\ Mw, M^o) 



a 

+ TT7 



-622(9', Mw, M^o) + 622(9', Mh+,Mao) (55) 



+ 4,622(9', Mh+,Mho) + c%J)22{q\ MH+,Mho) 
+ 4,622(9', Mw, Mho) + sij)22{q\ Mw, M^o) 



where the functions 60 and 622 are defined in eqs. (|66| , |67D , 4a — sin {(3 — a) 
and a is the mixing angle between the two MSSM scalar Higgs bosons 



defined e.g. as in ||31| , |30| . We recall that we use M^o, Mho and a with 
1-loop corrections included via the effective potential approach of ref. [^ . 
Each generation of sleptons contributes 



a 



47T{Uww{q'))'' = 2- cib22{q\ME,,M,) + 4622(9', M^j,, M^) 



(56) 



where cl = cos 0/, and the mass eigenstates of charged sleptons are defined 
by their relations to the left and right handed states: 



e; 



clE^ - slE^ 



E 



R 



slE^ +CLE2 



(57) 
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Similarily, one generation of squarks gives 



MUwwiq')) 



2\\sq 






44^22 (g', Mu„MdJ + 44^22 (g', Mu„Md,) 

+ 44^22(g', Mu„Md,) + 44^22(g', Mu„Mn,)] (58) 
with cu D = cos 0(7/5 defined by 



f/ 



Ct/f/f - St/f/2^ U^ = suUt+cuU^ 



Dl = cdD^ - sdD^ D^ = sdD-^ + cdD:^ (59) 

Finally contribution of two charginos and four neutralinos reads: 



47r(n^H^(g^)) 



2\\CN 



a 
s 



2 4 



Tj E E (I <c/iv I ' + I (^ic/N I ' ) [-4^22 (g' , mc^ , mjv, ) 



{<f - "^L - "^^J&o(g^ mc^,mNj 



(60) 



2 4 



^ E E 4^e(4c/Jvc4"c•/iv)"^c,"^7v.&o(g^ "^c, , "^jvJ 
^ j=ii=i 



The explicit formulae for the left and right couplings c'l rc/n '^^ charginos 
and neutralinos to W^ will be given shortly. 

For the tlzz the difference of the MSSM and SM Higgs sectors reads: 



2\\h 



47r{Uzz{q')) 






h{q\Mz,M^o) 



4,6o(g^ Mz, Mho) - sIMq^ Mz, M^o) 



+ 



a 



slJ)22{q\ Mz, Mho) + clJ)22{q\ Mz, Mho) (61) 



+ 4J)22{q^ Mao, Mho) + slj)22{q^ Mao, Mho) 



+ 



a 



p(? 



c' - rs')%2iq, Mh+,Mh+) - b22iq, Mz, M^o) 



One generation of sleptons yields: 



2\\sl 



AniUzziq')) 



a 



s^c^ 



b22{q\ M,, M,) + (4 - 2e)%22{q\ Me,, Me,) 

+ {sl-2s''fh22{q\ME„ME,) + 2clslh22{q\ME„ME,)\ (62) 
Likewise one generation of squarks contributes: 



47r(n^z(g')) 



2\\sq 



a 



Pc^ 



2cfjsijb22{q\ Mu„Mu,) + 244622(9', Md^Md,) 
+ (4 - fYb22{q\ Mu„ Mu,) + (4 - fs')%22{q\Mu„Mu;, 



+ (4 - -ffb22{q\Mn„Mn,) + (4 - -f)%22{q\MD„MD;, 



(63) 
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Two charginos and four neutralinos give 



47r(nzz(g')) 



2\\CN 



1 a 



2 2 



EE(l<cP + K,c\') -^h2{q\mc,,mc, 



'-2 2 



{q^ - ml - m^J6o(g^ mc,, mcj 



^J2J2 T^4c^L,c(^R*c)^c,mcM(l^^ ^c, , mc,) 



j=i i=i 



+ 



1 a 



4 4 



1 s2c2 



2 _2 _2 



- (g^ - "^?^, - "^^, ) &0 (g^ mN^ , mjv, ) 



(64) 



1 a 



4 4 



9^ EE T^^K,Nc'R*N)^N,mNMQ^ ^ mN^,mN,) 



2s 



j=i i=i 



The function 622 can be easily expressed in terms of standard a and 
60 functions [^, [l^ as: 



^22(g^"^l,"^2) 



-a(mi) + -a(m2) + -m^ bo (q"^, 7111,1712) 
b b 6 



2n2 



12^ 



,2_2 



-bo{q'^,'mi,m2) 
1 



(65) 



"^T-m2 . , . , i 2 , 2 1 2N 



-^2 2 [^("^1) ~ ^("^2)] + -{m^ + m2--q) 



Here a and 60 are the "renormahzed in MS scheme" functions: 

2 
a{m) =m^(-l + log-^) 



(66) 



7/2 N ; , 1 x(x - l)g^ + (1 - x)m? + xmn .^„- 

6o(g', mi, m2) = - / dx log ^ ^^^ '^2 ^ (67) 



In terms of the mixing matrices Z± and Zjq the couplings needed in 
eqs. (60,64) take the form I^O] |^: (the fermion i is incoming and j is 
outgoing) 



ci 



L,C/N 



c' 



R,C/N 



ylj* y2i -*- r72j* ryii 



(68) 



Factors cj^'^ and c^ ^ given in BQ] are not properly symmetrized. 
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'^L,C — 


z^^'z^; + [& - s^w' 


Ji — 
^R,C — 


Z^jZ^_!* + (c^ - s^)6^' 






Ji 





(69) 

and 



(70) 

We end this Appendix with two comments. First, as long as neutrahnos 
appear only in closed loops it is not necessary to do anything with negative 
eigenvalues of their mass matrix (53). Using {^ + m) / {p"^ — m"^) for m<0 
as a propagator gives the same physical results as other more complicated 
prescriptions which require modification of couplings (see eg. 0]). It is 
possible, however, to choose the Z^- matrix such that all mass eigenvalues 
are positive. 

Second, the transition from Weyl spinors Aj ' to four component fields 
ipN- and tp(j+ is not unique. One can work for example with 



^c- = T+ (71) 




instead of 1^^+- This freedom in building four component fields from two 
component ones is in fact the basis of the recently developped technique ^3 



for handling Majorana particles which often produce "clashing arrows" in 
Feynman diagrams. However realizing that (at the level of lagrangian) e.g. 

a^%Al - 7')^^.0+ = a^V^,(l - 7')%0+ (72) 

and that in different vertices of a given Feynman diagram different forms of 
the same interaction can be used allow avoiding essentially all problems with 
Majorana particles. 
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FIGURE CAPTIONS 

Figure 1. 

Comparison of sin^ e{Mz)FULL with sin'^ e{Mz)LLT in the MSSM for 
sparticle spectra corresponding to the general scan described in the text in 
four different cases: 

a) mt = 160 GeV tan/? = 2, 

b) mt = 160 GeV tan/? = 50, 

c) nit = 180 GeV tan/3 = 2, 

d) mt = 180 GeV tan/? = 50. 

Cirles denote spectra with all SUSY particles heavier than 500 GeV and 
50 < Mao < 500 GeV. Stars correspond to spectra with M^o < 250 
GeV, rriCi, "^a^ — 250 GeV and heavy sfermions, > 500 GeV. Squares 
correspond to spectra with all sparticles light, < 250 GeV and 50 < M^o < 
250 GeV. 

Figure 2. 

a) sin^ 9{Mz) and b) a{Mz)EM in the MSSM as a function of the (com- 
mon for all the three generations) sneutrino mass. All other SUSY particles 
(including the A^ Higgs boson) have masses 2 TeV, with no Left - Right 
mixing. Solid lines correspond to the full calculation, dashed ones to the 
LLT approximation, dot-dashed lines show our approximation (eqs. ( [HT|P^ ) 
combined with the fit (pS])). Two dotted lines show the results of application 



of the Faraggi Grinstein formulae (see |Ty], eqs. (5.10)) to the SM couplings 



obtained from our fit (p5|) 



Figure 3. 

a) sin^ 6'(M^) and b) a{Mz)EM in the MSSM as a function of the (purely) 
left handed sbottom mass. All other SUSY particles (including the A^ Higgs 
boson) have masses 2 TeV, with no Left - Right mixing. Solid lines corre- 
spond to the full calculation, dashed ones to the LLT approximation, dot- 
dashed lines show our approximation (eqs.(^,^) combined with the fit (pS])). 
Two dotted lines show the results of application of the Faraggi Grinstein for- 
mulae (see [jlO[, eqs. (5.10)) to the SM couplings obtained from our fit (|25|). 



Figure 4. 

a) sin^ 6'(M^) and b) a{Mz)EM in the MSSM as a function of the pseu- 
doscalar mass M^o. mass. All SUSY particles have masses 2 TeV, with no 
Left - Right mixing. Solid lines correspond to the full calculation, dashed 
ones to the LLT approximation, dot-dashed lines show our approximation 
(eqs.(0,^) combined with the fit (^5]) ). Three dotted lines show the results 
of apphcation of the Faraggi Grinstein formulae (see |TD[, eqs. (5.10)) to the 
SM couplings obtained from our fit (pSj). 
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Figure 5. 

Parameter /x as a function of Tsusy for spectra of ref. p4|, |25[ obtained for 
different values of nit and tan/? and subject to the conditions: i) Squark 
masses below 2 TeV, ii) Sparticle masses above the experimental limits. 

Figure 6. 

Values of a^ and s^ (vertical and horizontal axis respectively) calculated 
with the use of sparticle spectra characterized in Fig. 5 and for the same rrit 
and tan/? values. Stars (diamonds) denote LLT (full) calculation of the 
MSSM couplings. The upper (lower) group of points in Figs, a), c) and d) 
corresponds to af^^ = 0.125(0.115). The three groups of points in Fig. 
b) correspond to a^'^^ = 0.13, 0.125, 0.115 .Contours in the plot denote 
bounds imposed by SUSY SU(5) unification conditions with dem (in the 
MSSM) allowed to vary in the range: 1/129.0 - 1/131.8: Solid lines - no 
GUT scale thresholds; dashed (dash-dotted) lines - M^^^/M™*^ = 30 (1000) 
(for details see text). The three dotted hues correspond to fixed values of 
Mhs = 1 X 10^^ (uppermost), 2 x 10^^ 5 x 10^^ (for Aem = 1/131.8). The 
region below the long-dashed lines correspond to violation of perturbativity 
condition Fj2(MGc/r)/47r < 1. 

Figure 7. 

Values of af^ predicted by SUSY unification without GUT scale thresh- 
olds as a function of Tsusy calculated with the use of sparticle spectra 
characterized in Fig.5 and for the same rrit and tan/? values. Stars 
(diamonds) denote LLT (full) calculation of the MSSM couplings. The re- 
gion below the long-dashed line corresponds to violation of perturbativity 
condition (Yj2(Mg[/t) /47r < 1). 

Figure 8. 

Mh3 as a function of Tsusy for df^^(Mz) = 0.125 and sparticle spectra 
characterized in Fig.5. Stars (diamonds) denote LLT (full) calculation of the 
MSSM couplings. 
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